Rydberg transition frequencies from the Local Density Approximation 
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A method is given that extracts accurate Rydberg excitations from LDA density functional cal- 
culations, despite the short-ranged potential. For the case of He and Ne, the asymptotic quantum 
defects predicted by LDA are in less than 5% error, yielding transition frequency errors of less than 
O.leV. 
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The excitation energy spectrum of atoms, molecules, 
clusters, and solids can now be accurately calculated via 
Time-dependent Density Functional Theory (TDDFT) 
0. In Casida's matrix formulation 0, first the self- 
consistent solution of the ground-state Kohn-Sham (KS) 
equations is found Q, and the differences between oc- 
cupied and unoccupied KS orbital energies may then be 
regarded as a first approximation to the true excitations 
of the system. In a second step, these KS frequencies are 
corrected to become the true transitions of the many- 
body system. The quality of the results depends cru- 
cially on the functional employed for the solution of the 
ground-state problem. 

The Local Density Approximation (LDA) is the sim- 
plest and historically most successful approximation in 
DFT 3]. Whereas new generations of functional have 
achieved better accuracy than LDA for many properties, 
its ratio of reliability to simplicity has no paragon. But a 
well-known shortcoming of the ground-state LDA poten- 
tial of an atom or molecule, already recognized by Tong 
and Sham in the early days of the LDA |4j], is that it 
decays exponentially at large distances, rather than as 
— 1/r as the exact KS potential does. As a consequence, 
the LDA potential does not support a Rydberg series of 
bound states. Also, the magnitude of the highest occu- 
pied molecular orbital (HOMO) is typically too small by 
several eV in an LDA calculation, so HOMO — > Ryd- 
berg transitions appear as HOMO — > continuum excita- 
tions. Because of this limitation, several schemes have 
been devised to asymptotically correct the LDA poten- 
tials. Casida and Salahub (CS) [f| proposed a "shift-and- 
splice" approach consisting on shifting the LDA potential 
downwards in the bulk regions and joining it continuously 
with the van Leeuwen-Baerends potential [(| where they 
cross. Tozer and Handy (TH) proposed a very similar 
procedure using a Fermi- Amaldi tail as suggested by 
Zhao, Morrison and Parr and a more sophisticated 
way to smoothen the transition from the inner to the 
outer potential. Wu, Ayers, and Yang (WAY) Q, not- 
ing that neither of these potentials yield an energy min- 
imum, proposed a variational method for correcting the 
exchange-correlation (xc) potential. Their construction 
imposes the correct asymptotic behavior by using the 
Fermi-Amaldi form as a fixed reference potential, and 
the coefficients of a linear combination of basis functions 



are determined through minimization of the energy for 
a given choice of energy functional. Their asymptoti- 
cally corrected LDA potentials very closely resemble the 
pure LDA ones shifted downwards in energetically im- 
portant regions. For any finite basis set, however, they 
are no longer functional derivatives of the LDA energy 
functional. Interestingly, the results obtained so far for 
excitation energies with the WAY method are not as good 
as those obtained with TH potentials that retain the pure 
LDA form in the core regions ^J. All of these methods 
have improved upon LDA on the prediction of Rydberg 
excitation frequencies, leading to the conclusion that the 
correct long-range behavior of the potentials is indispens- 
able to describe such excitations. But we show here that 
in fact it is not. 

The purpose of the present work is to show how the 
Rydberg excitation energies are encoded in the short- 
ranged LDA potentials. Imposing the correct asymptotic 
behavior is one way to decode them, but not the best 
one, since different tails lead to different answers (e.g. 
van Leeuwen-Baerends vs. Fermi-Amaldi) , obscuring the 
information provided by a pure LDA calculation. 

It was recently shown [llj that in spite of incorrectly 
describing photoabsorption as if it were photoionization, 
the oscillator strengths of Rydberg excitations show up 
in the LDA spectrum as continuum contributions with 
excellent optical intensity. The dipole matrix elements 
for HOMO — > Rydberg transitions are accurate in LDA 
because (1) the shape of the LDA HOMO is very close to 
that of the exact KS HOMO, even if its energy is not, and 
(2) LDA continuum orbitals at frequencies corresponding 
to HOMO — > Rydberg transitions, are also very close to 
the exact KS Rydberg orbitals in the crucial region for 
optical absorption, i.e., where the HOMO has high am- 
plitude. The underlying physical reason for this is simple 
[llj: the LDA xc-potential is very close to the exact xc- 
potential near the nuclei, and runs almost parallel to it 
in the valence regions |12j . 

An obvious objection to this claim of success of the 
LDA is that it cannot predict the positions of the Ry- 
dberg excitations, even if it produces an ionization en- 
velope that approximates well the discrete photoabsorp- 
tion spectrum. We now show that, in fact, LDA does 
predict the position of high-rt Rydberg excitations very 
accurately. We use concepts of quantum defect theory, 
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developed before the advent of DFT by Ham [13| and 
Seaton [2(|. The quantum defect n n \ parametrizes the 
energy E n i of a Rydberg state as: 
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For an electron orbiting in a Coulomb field outside an 
ionic core, as in a high-rt Rydberg state, fj, n i represents 
the effect of the field that prevails within the core. Al- 
though the Coulomb field outside the core is invoked for 
its definition in Eq.Q, the actual number /i„; is deter- 
mined only by the forces within [lfij . It is typically a 
very smooth function of n, and approaches rapidly the 
asymptotic quantum defect, [tool, as n — > oo. 

We will focus on the KS asymptotic quantum defect 
of the (I — 0)-Rydberg series that converges to the first 
ionization threshold of an atom (the I = subscript will 
be dropped from now on). We first propose a method 
to extract \i n from a given orbital, and illustrate it with 
a simple example. This method was inspired by Fano's 
original discussion in Ref . [Til] . We then apply it to the 
cases of He and Ne, where the exact KS quantum defects 
are known, and show that the LDA produces /Zoo's which 
are in less than 5% error. 

Quantum defect from an orbital: Consider a long- 
range potential that equals —1/r for r > ro. The so- 
lution of the radial Schrodinger equation is well known 
for r > tq. For negative energies E < 0, the physically 
acceptable solutions are Whittaker functions (we will re- 
strict the analysis to s-states): 



(f> >ro (r) = AW 1/kA/2 (2kr) 



(2) 
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FIG. 1: s-quantum defect as a function of n for the potential 
shown in the inset. The quantum defect is a smooth function 
of n and converges rapidly to its asymptotic value, fi^ = 
-0.441. 




quantum defect = -0.441 



FIG. 2: Solution of Eq.© for /i as a function of r, for the 
potential of Fig0 the logarithmic derivative of the n = 20 
orbital was found numerically as a function of r, and at each 
location r, it was inserted into Eq.lJ^J to get fi(r). 



where A is a constant and k = y/2\E\. The logarithmic 
derivative of 4>>r is given by: 

d\n(j> >ro _ I n* 1 U(-n*;2;2r/n*) 

dr n* r rU(l-n*;2;2r/n*) ^ ' 

Here k was written as k — (n*) _1 , with n* — (n — fx n ), 
where n numbers the bound state, and [i n is the quantum 
defect; U is the confluent hypergeometric function |1(t| |. 

Regardless of the shape of the potential for r < tq, the 
logarithmic derivative of 4>< ro must equal that of 4>>r at 
ro. Now suppose that an orbital is given to us, with the 
information that it is the n — 8 state of a potential that 
possesses a Coulomb tail. We can immediately obtain 
H from this orbital by solving Eq.Q numerically, using 
n = 8 and some large value of r. If we observe that n(r) 
changes as r is increased, we can conclude that Eq.@ is 
being used in the region where r < fo, and its solution 
cannot be interpreted as the quantum defect. 

For example, consider a potential which is equal to a 
constant C for r < ro and to —1/r for r > tq. For ro = 1 
the matching condition is: 

r , r 1 * U{-n*;2;2/n*) 

k coth k = n* - —i ' ' . '-. 4 

n* £7(l-n*;2;2/n*) v ' 



where k = \J2\E - C\. Eq. |@J| was solved for the first 20 
bound states. For C = ro = 1 the asymptotic quantum 
defect is [i^ = —0.441, see FigQ] Figure [21 shows /«(r), 
the solution of Eq.(j2J as a function of r, for the n = 20 
orbital. Clearly, the quantum defect for a given state can 
be obtained by looking at the respective orbital anywhere 
in the region r > rg. In particular, it can be obtained at 
ro. It represents the accumulation of phase due to the 
non-Coulombic potential in the region of r < ro up to 
r = r . 

Imagine now that the potential is altered by truncating 
the Coulomb tail far away, at T\ » ro, and making 
the potential equal to the constant — l/r\ for all r > 
n. This modified potential has an orbital which -up to 
a constant- is almost identical to the original n = 20 
Rydberg orbital in the region r < ro, but it may now 
be a scattering orbital, behaving very differently in the 
region r > t\ . We can still solve Eq. J^) on this scattering 
orbital, with n = 20 on the right-hand side, and find 
fi ~ —0.441 at ro. The altered potential does not have 
a Rydberg series, yet the solution of Eq.@ at r can 
still be interpreted as the asymptotic quantum defect of 
the Rydberg series that was lost as a consequence of the 
alteration. 
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FIG. 3: s-quantum defect as a function of n for the exact KS 
potential 11711 of the He atom. The quantum defect converges 
rapidly to its asymptotic value, /loo = 0.213. 
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The Coulomb tail has nothing to do with the value of 
fi. It only has to do with its definition. 

One technical point deserves comment: rather than 
giving /i(r) directly, Eq.© yields n*(r), and it would 
seem that it needs to be solved with ever increasing ac- 
curacy as n — > oo, since orbitals of different n become es- 
sentially identical at large n in the [0, ro] region. Only the 
fractional part of the solution of Eq.@ is to be trusted, 
yielding the fractional part of n* , which is also equal to 
the fractional part of the asymptotic quantum defect, de- 
noted {fJ.oo}- The integer part of fXoo can be determined 
by a simple node counting as \pLoo] — N — Nc, where N is 
the number of nodes of the given orbital in the [0, ro] in- 
terval and Nc that of the correspoding pure Coulomb or- 
bital. The asymptotic quantum defect /ioo = [/i o]+{Atoo} 
is thus fully determined this way. 

Results for Helium and Neon: Consider the first s- 
Rydberg series of the He atom. Figure shows the s- 
quantum defect as a function of n from the exact KS 
potential obtained by Umrigar and Gonze 17]. The 
asymptotic quantum defect /Zoo = 0.213 can be extracted 
through Eq.© from e.g. the n = 20 orbital, just as it 
was done in the previous example. It is clear from the 
solid line of FigQ]that at ro ~ 1, Eq.JJJ is already giving 
an accurate value of /i, which is quite remarkable consid- 
ering that the KS potential at r ~ 1 is still not equal to 
— 1/r (see FigEJ . Most of the quantum defect is built up 
close to the nucleus (steep rise for < r < 1 in FigQ}, 
and its final value has been reached before the potential 
becomes purely coulombic. 

The LDA potential runs almost parallel to the ex- 
act one in the region 1 < r < 2 (where /i m can al- 
ready be extracted accurately), and orbitals correspond- 
ing to the same frequency (exact and LDA) are there- 
fore very close in that region, see FigHj] In the spirit 
of Ref . Jljj , we compare the exact energy- normalized 20s 
orbital (which is essentially identical to the zero-energy 
state in the region < r < 6) and the LDA orbital of en- 
ergy I + e^ s DA = 0.904 - 0.571 = 0.333. Notice how good 
the LDA orbital is in the region 1 < r < 2. We show in 
FigQJthe solution of Eq.J3J) when this scattering LDA or- 
bital is employed. Clearly, the plateau of the LDA curve 
in the 1 < r < 2 region is an accurate estimate of the 
quantum defect. The value of /i on this plateau is 0.205, 
an underestimation of less than 4% with respect to the 



FIG. 4: He atom: solution of Eq.© for /i as a function of 
r; The n = 20 orbital was used for the exact case, and the 
scattering orbital or energy E — I + €i^ A was used for the 
LDA. 




-2 -1.5 -1 -0.5 0.5 1 
ln(r) 



FIG. 5: Comparison of the exact KS potential of the He atom 
[l^l (solid line), and the LDA potential (dashed line). The 
Coulomb potential is also shown (dotted line). At r ~ 1 the 
exact potential is almost Coulombic. The inset shows the 
potentials themselves. 



exact value. 

Thus, given the ionization potential of the system, 
LDA gives a very accurate prediction of the asymptotic 
quantum defect. The ionization potential is needed to 
choose the appropriate LDA scattering orbital, but the 
results are not terribly sensitive to it. We repeated the 
same procedure with the LDA ionization potential (de- 
fined as £ L DA(He)-£: LD A(He + )=0.974) instead of the 
exact one, and found /J^ A — 0.216, overestimating the 
exact by just 1%. 

Our analysis provides a natural way to asymptotically 
correct the LDA potential: simply force the LDA plateau 
of FigQ] to stay constant for all r. The resulting func- 
tion n(r) determines a zero-energy orbital which in turn 
uniquely determines a potential through inversion of the 
Kohn-Sham equations at zero energy. We emphasize, 
however, that such potential is not needed to obtain /z^, 
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FIG. 6: Radial orbitals of He: LDA orbital of energy E=0.333 
and exact 20s orbital; the HOMO is also shown. The LDA 
scattering orbital and the exact KS Rydberg orbital are very 
close in the region 1 < r < 2, where the quantum defect can be 
extracted (see text, and Fig^J. The fact that the LDA orbital 
has an incorrect asymptotic behavior at large r is irrelevant 
for the value of fi, as well as for optical absorption |Tl|. 



TABLE I: Transition frequencies (in eV) for the first six 
discrete 2p — » ns transitions in Ne, from the exact and LDA 
KS potentials. 



trans. 


transition frequency 
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exact 
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3s 


16.468 


16.604 


2p- 


4s 


19.603 


19.666 
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20.534 


20.561 
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6s 


20.931 


20.945 
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7s 


21.135 


21.143 


2p- 


8s 


21.255 


21.260 



but the question of what long-range potential gives rise 
to the same /ioo is certainly worth addressing. The TH Q 
or CS @ methods for asymptotically correcting the LDA 
potentials require choosing a radius ro where the proper 
tail is to be pasted. Our analysis also provides a way to 
rigorously justify such choice, since rp should clearly be 



on the LDA plateau of /x(r). Minimizing d^i(r)/dr\ ro de- 
termines its precise value of r ~ 1.3. To test this, we per- 
formed a simple Latter-type asymptotic correction 
by pasting a—l/r tail to the LDA potential shifted down- 
wards by Ulda(^o) ~ r o : an d scanning through ro. The 
errors in both fi^ and e\ s were minimized at ro ~ 1.5. 

Repeating the same procedure for the Ne atom we 
found again a distinctive plateau in the LDA curve of 
^oo(r) at ^J^P A = 1.366, an overestimation of 4% with 
respect to the exact value (/ti^ act = 1.313). 

It has been shown [l^ for the case of Ne that 3 num- 
bers are enough to fit very accurately the entire curve 
of Kohn-Sham quantum defects of a given /. Although 
we have only determined one such number here (/^oo), we 
use it nonetheless to approximate all the / = KS or- 
bital energies as e^nA™ = ~i^( n ~ Moo)] -1 ' The results 
are presented in Table JT| The errors are smaller than 
typical TDDFT errors jj. 

Finally, our results suggest that the LDA can also be 
employed to calculate accurate low-energy electron-ion 
scattering phase shifts. In fact, the LDA asymptotic 
quantum defect found in this work, /x^ A = 0.205, im- 
mediately yields, through Seaton's theorem J2jj, a pre- 
diction for the zero-energy s-phase-shift for electron-Hc + 
scattering: S(E = 0) L da = = 0.644. The exact 

KS phase shift is 5(E = 0) CX act = 0.669. This value is 
also remarkably close to the average of the experimental 
singlet/triplet zero-energy phase shifts |2ll |. 

In conclusion, we have shown that rather than mod- 
ifying the shape of the LDA potentials one can modify 
the interpretation of the results of a pure DFT-LDA cal- 
culation. The results are excellent for the excitation to 
high-lying s-Rydberg states in He and Ne, and we are 
working to extend these ideas to Rydberg excitations in 
molecules. 

This work was supported by NSF grant No.CHE- 
0355405. 
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